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We show that measurements of the forward-backward charge asymmetry {Afb{M, y)) of Drell- 
Yan dilepton events produced at hadron colliders provide a new powerful tool to constrain Parton 
Distribution Functions (PDFs). PDF uncertainties are the dominant source of systematic error in 
precision measurements at hadron colliders (e.g. sin^ 9eff{Mz), sin^ 9w = I and the 

mass of the W boson). We show that the values of fits to extract sin^ 9l'jj{Mz) from A/rB(M,y) 
with different PDF replicas can be used to place additional constraints on PDFs. In turn, using 
these constrained PDFs significantly reduces the PDF errors in precision measurements of elec- 
troweak parameters. The measurement of the on-shell sin^ 9\y = I /M| is equivalent to an 
indirect measurement of the W mass. The errors in this indirect measurement of the W mass are 
competitive with direct measurements. For example, with 200 fb^* at 13 TeV, the expected error 
in the indirect measurement of the W mass is ±9 MeV. 
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Within the standard model, measurements of the mass of the Z boson and top quark, in com¬ 
bination with the mass of the Higgs boson, can be used to predict the mass of the W boson. At 
present, the average of all direct measurements of the mass of the W boson (80385±15 MeV) is 
about one standard deviation higher than the prediction of the standard model. Predictions of su¬ 
persymmetric models for the W mass are also higher than the predictions of the standard model. 
Therefore, more precise measurements of the mass of the W boson are of great interest. 

Alternatively, the W mass can also be extracted indirectly from measurements of the on-shell 
electroweak mixing angle sin^ dw by using the relation sin^ 6vv = 1 ~ jM^- Measurements of 
the forward-backward charge asymmetry in Drell-Yan dilepton events produced at hadron colliders 
(in the region of the Z pole) have been used to measure the value of the effective electroweak 
(EW) mixing angle sin^ 0^jy^(Mz)[l, 2, 3, 4]. In addition, by incorporating electroweak radiative 
corrections in the analysis the CDF collaboration has also measured the on-shell EW mixing angle 
sin^6vv[h 2]. An error of ±0.00030 in the measurement of sin^Ovy is equivalent to an indirect 
measurement of the W mass to a precision of ±15 MeV. However, the PDF error quoted in the 
most recent measurement of sin^ deff by the ATLAS collaboration [4] at the LHC is ±0.00090. 
Therefore, a significant reduction in the PDF errors is needed. Here, we show how Afb data (both 
at the Tevatron and LHC) also provide a new powerful tool to constrain PDFs. Addition details can 
be found in Ref. [5]. 





Figure 1: Left panel: Afb{M) at the Tevatron for uu (black), dd (red) and the sum of the two (blue). Right 
panel: (a) The difference between Afg(M) at the Tevatron for 10 NNPDF3.0 replicas and AfB{M) calculated 
for the default NNPDF3.0. Here sin^ 9w is fixed at a value of 0.2244. (b) The difference between Afb{M) 
for different values of sin^ Ow ranging from 0.2220 (shown at the top in red) to 0.2265 (shown on the bottom 
in blue), and Afb{M) for sin^ 0vr=O.2244. Here, the NNPDF3.0 default PDF is used. 


For pp collisions, the direction of the quark is predominately in the proton direction, and the 
direction of the antiquark is predominately in the antiproton direction. Here, most of the cross 
section originates from the annihilation of quarks in the proton with antiquarks in the antiproton. 
Therefore, Afb is measured under the assumption that the quarks originate form the proton, and 
the antiquarks originate from the antiproton. 
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The extraction of sin^ from Afb{M) at the Tevatron is sensitive to PDFs for two reasons. 
First, Afb{M) for up and down type quarks is different as shown in Fig. 1. The asymmetry at 
the Tevatron originates primarily from up quarks and is diluted by the fraction of down quarks in 
the proton because the asymmetry for down quarks is much smaller. In addition, there is a small 
fraction of events for which the annihilation is between sea antiquarks in the proton with a sea 
quarks in the antiproton. This also results in a dilution of the measured asymmetry. 

The mass dependence of Afb{A1) depends on both sin^ dw and on PDFs In the region of the 
Z pole, Afb{M) is sensitive to the vector coupling, which depend on sin^ dw- At higher and lower 
mass AfB{M) is sensitive to the axial coupling and therefore insensitive to value of sin^Pvv- In 
contrast, the magnitude of the dilution of Afb{M) depends on the PDFs. The sensitivity to PDFs is 
largest in regions where Afb{M) is large (i.e. away from the Z pole). 
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Figure 2: LHC; Left-Top panel- Afb at the LHC at ^/s=S TeV for six rapidity bins. The horizontal scale for 
each of the six plots is the invariant mass. Left-Bottom panel; The green bands span the difference be¬ 
tween A/7B(M,y) calculated for the 100 NNPDF3.0 replicas and calculated for the central default 

NNPDF3.0 for the six rapidity bins. The blue lines are the differences between AFB{M,y) calcu¬ 

lated with different values of sin^ 9gff and the values calculated with nominal sin^ 0gy^^^=O.2312O. Right-Top 
panel. Analysis of one of the 64 LHC pseudo experiments The top two panels show the extracted sin^ Be/f 
and corresponding Xafb values from fits to AFB{M,y) versus replica number for the 100 NNPDF3.0 replicas. 
The Right-Bottom panel shows the same results in the form of a scatter plot of Xafb values versus sin^ Be// 
for one pseudo experiment. 

The right panel of Fig. 1 shows the sensitivity of Afb{M) at the Tevatron to PDFs. Also 
shown is the sensitivity of Afb{AI) at the Tevatron to sin^Fvv- There is a large difference in the 
Afb{AI) predictions for PDF sets with different ^(x) and ^{x) in regions where Afb{M) is large 
and positive (M>100 GeV). The changes in Af 5 (M) in regions where A/7 b(M) is large and negative 
(M<80 GeV) are the opposite direction. In contrast, different values of sin^ dw change Afb{M) 
primarily in the region near the Z pole. However, here the change is in the same direction above 
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and below the Z pole. Therefore, if we extract sin^ dw from AfB{M) data using different PDFs, 
PDFs with poor values of are less likely to be correct. 

The Left-Top panel of Fig. 2 shows AfB{M,y) at the LHC at v^=8 TeV for six rapidity bins 
(i=l to 6) with average y values of 0.2, 0.6, 1.0, 1.4, 1.8 and 2.2. The horizontal scale for each of the 
six plots is the invariant mass. The calculations are done with the POWHEG MC generator. 

The version of POWHEG that is used does not include electroweak radiative corrections. There- 

'y 

fore, POWHEG requires an input value of sin for the calculation of Afs- The green bands 
span the difference between AfB{M,y) calculated with the 10 NNPDE3.0 replicas and AfB{M,y) 
calculated with the default NNPDE3.0 PDE. The blue lines are the differences between AfB{M,y) 
calculated for several values of sin^ and Afb{^a) for th® nominal sin^ deff=f).23\20. Eor all 
of the blue lines, AFB{M,y) is calculated with the default NNPDE3.0 PDE. 

At the EHC, as for the Tevatron, the dependence of AfB{M,y) on sin^de/f and on PDEs is 
different. In the region of the Z pole, AfB{M,y) is sensitive to the vector coupling, which are 
functions of sin^ %/• At higher and lower mass AfB{M,y) is sensitive to the axial coupling and 
therefore insensitive to value of sin^ As is the case for the Tevatron, the magnitude of the 
dilution of AfB{M) is larger in regions where the absolute value of Afb{A1) is large (i.e. away from 
the Z pole). At the EHC the dilution depends on both M and y. The combined mass and rapidity 
dependence of the dilution at the EHC provides more stringent constraints on PDEs than Afb{M) 
measurements at the Tevatron. 

The NNPDE3.0 PDE set is given in the form of N (e.g. 100 or 1000) replica PDEs. Each of 
the PDE replicas has equal probability of being correct. The central value of any observable is the 
average of the values extracted using each one of the N PDE replicas. The PDE error is the RMS 
of the values extracted using each of the N replicas. 

One advantage of the PDE replica method is that constraints from new data can easily be 
incorporated in any analysis by using different weights for each replica. Replicas for which the 
theory predictions are in agreement with the new data are given higher weights, and replicas for 
which the predictions are in poor agreement are given lower weights. The weights are derived from 
the x^ values of the comparison between the new data and theory prediction using each of the PDE 
replicas. The central value of any observable is then the weighted average of the values extracted 
using each one of the N PDE replicas. The PDE error is the weighted RMS of the values extracted 
using each of the N replicas. 

The procedure of constraining a PDE set with new data was initially proposed by Giele and 
Keller[8]. They proposed that each of the N PDE replicas be weighted by w;, and the weights 
reduce the effective number of replicas[l 1] from N to Neff - Here 

e-iXf 1 N 

More recent discussions of the method can be found in references [6, 7, 9, 10, 11]. The mass and 
rapidity dependence of Afb can be used to both provide additional constraints on PDEs and reduce 
the PDE error in measurements of sin^ dw- 

Eor studies of AFB{M,y) at the EHC we simulate Drell-Yan data for 64 pseudo exper¬ 

iments for a CMS like detector at v^=8 TeV. The pseudo data is generated using the POWHEG 
NEO MC generator with the default NNPDE3.0 PDEs and sin^ =0.23120. 
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Table 1: Values of sin^ 9w with statistical errors and PDF errors expected for a 15 fb^* Drell-Yan 
sample at the LHC (at 8 TeV). The pseudo data is generated by the POWHEG MC generator with the default 
NNPDF3.0 PDF, and sin^ 0g^y^=O.2312O. The PDF error for a standard analysis is compared to the PDF 
error for an analysis with both x\fb weighting and Xafb + ^Wasym weighting. In addition, expected errors 
for larger statistical samples are shown. 


input 

POWEG 

Default 

NNPDE3.0 

(261000) 

EHC CMS like 
Pseudo-Exp. 

15 fb-' 8 TeV 

6.1M 

reconst, events 

EHC CMS like 
Pseudo-Exp. 

19 fb-' 8 TeV 
15M e^e^) 

reconst, events 

EHC CMS like 
Pseudo-Exp. 
200fb-i 13 TeV 

120M 

econst. events 

sin^ deff statistical error 
sin^ deff CTIO PDE error 

±0.00050 

± 0.00080 

±0.00034 

±0.00011 

NNPDE3.0 Average 
PDE error RMS 

Neff=m 

±0.00051 



xIfb weighting 
weighted PDE error RMS 

Neff = 46 

±0.00029 



XAFB+Xwsym Weighting 
weighted PDE error RMS 

Neff = 21 

±0.00026 

± 0.00022 

±0.00014 

Asin^ deff Stat-i-PDE 
AMw indirect Stat-i-PDE 

±0.00056 

±28 MeV 

±0.00040 

± 20 MeV 

±0.00018 

±9 MeV 


For each pseudo experiment, we generate a sample of 15.6 Million ^ events with > 
50 GeV, which corresponds to an integrated luminosity of 15.0 This is similar to the k,\9 
fb^' of integrated luminosity collected by CMS and ATLAS at 8 TeV. To this sample, we apply 
acceptance and transverse momentum cuts which are similar to a CMS-like detector. We also 
smear the events with a muon momentum resolution similar to a CMS-like detector. The final 
sample consists 6.7M reconstructed events. 

The 8 TeV W asymmetry data at the LHC has not yet been incorporated into the most recent 
PDF fits. Therefore, in addition to AfB(M,y), we also use the default NNPDF3.0 PDF to generate 
pseudo data for the W decay muon asymmetry as a function of muon rapidity (for muon transverse 
momentum PT>25 GeV). This simulates the W asymmetry measurement at 8 TeV. 

In the analysis of each of the 64 pseudo experiments generated with the default NNPDF3.0 
PDF the simulated values of Afe(M,y) for each experiment are compared to Afg(M,y) templates. 

r\ 

The templates are generated with the POWHEG MC for a range of values of sin Qeff for each 

'y 

of the 100 NNPDF3.0 PDF replica. For each replica we extract the best fit value of sin deff, the 
corresponding statistical error and the fit Xafb- 

In addition, we calculate Xwasym which is the x^ for the agreement between the predictions for 
the W lepton decay asymmetry and the W decay lepton asymmetry pseudo data at 8 TeV for each 
of the 100 PDF replicas. 

The right panels of Fig. 2 shows the results from one of the 64 pseudo experiments at the 
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LHC. The top two panels on the right show the extracted sin^ dgff and corresponding Xafb values 
from fits to ApsiM^y) versus replica number for the 100 NNPDF3.0 replicas. The bottom panel on 
the right shows the same results in the form of a scatter plot of Zafb values versus sin^ dgff for one 
pseudo experiment. For each pseudo experiment we find fhe mean value and PDF error of sin^ dgff 
from fhe average and RMS of fhe sin^ exfracfed values using each of fhe 100 PDF replicas. 
The average and RMS values are done in fhree ways: (1) Using fhe sfandard average and RMS of 
fhe sin^ fil values. This analysis resulfs in a sfandard PDF error of ±0.00051 wifh 100 replicas. 

(2) Using fhe Xafb values of fhe fifs fo AfB{M,y) fo form a weighted average and weighted RMS 
of fhe sin^ values. This analysis resulfs in a PDF error of ±0.00029 wifh 46 effecfive replicas. 

(3) Using fhe combined XAFB'^Xwasym the fifs fo Drell-Yan AfB{M,y) pseudo dafa and fhe fifs 
fo fhe W lepfon decay asymmefry pseudo dafa fo form fhe weighted average and weighted RMS of 
fhe sin^ values. This analysis resulfs in a PDF error of ±0.00026 wifh 21 effecfive replicas. 

Table 1 shows values of sin^ dw wifh sfafisfical errors and PDF errors expecfed for a 15 fb^t 
Drell-Yan sample af fhe LHC (af 8 TeV). The pseudo data is generated by the POWHEG MC 
generator with the default NNPDF3.0 PDF, and sin^ 0eyy=O.2312O. The PDF error for a standard 
analysis is compared to the PDF error for an analysis with both Xafb weighting and Xafb + Xwasym 
weighting. As shown in Table 1, the number of effective PDF replicas is reduced when we apply 
constraints from Xafb ^^tl Xwasym- Therefore, the analysis will be more robust if we start with 1000 
PDF replicas. 

Also shown are the expected errors for larger statistical samples. With larger statistical sam¬ 
ples, the PDF constraints are more stringent, and the PDF errors are also reduced. The errors in 
this indirect measurement of the W mass are competitive with direct measurements. For example, 
with 200 fb^^ at 13 TeV, the expected error in the indirect measurement of the W mass is ±9 MeV. 
Additional details and studies for both the Tevatron and LHC are give in ref. [5] . 
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